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Abstract. RX J1826.2-1450/LS 5039 has been recently pro- 
posed to be a radio emitting high mass X-ray binary. In this 
paper, we present an analysis of its X-ray timing and spec- 
troscopic properties using different instruments on board the 
RXTE satellite. The timing analysis indicates the absence of 
pulsed or periodic emission on time scales of 0.02-2000 s and 
2-200 d, respectively. The source spectrum is well represented 
by a power-law model, plus a Gaussian component describing a 
strong iron line at 6.6 keV. Significant emission is seen up to 30 
keV, and no exponential cut-off at high energy is required. We 
also study the radio properties of the system according to the 
GBI-NASA Monitoring Program. RX J1826.2-1450/LS 5039 
continues to display moderate radio variability with a clearly 
non-thermal spectral index. No strong radio outbursts have 
been detected after several months. 

Key words: Stars: individual: RX J1826.2-1450, LS 5039, 
NVSS J 1826 14- 145054 - X-rays: stars - Stars: variables: 
other - Radio continuum: stars 



1. Introduction 

The star LS 5039 is the most likely optical counterpart to the 
X-ray sourceRX J 1826.2— 1450. Such an association was orig- 
inally proposed by Motch et al. (1997), hereafter M97, as a 
result of a systematic cross-coiTelation between the ROSAT 
All Sky Survey (Voges et al. 1996) and several OB star cata- 
logues in the SIMBAD database. The unabsorbed X-ray lumi- 
nosity, at an estimated distance of 3. 1 kpc, amounts to ix(0. 1- 
2.4 keV) ~ 8.1 X lO''^ erg s^^, and the hardness of the source 
is well consistent with a neutron star or a black hole, accret- 
ing directly from the companion's wind (M97). In the optical, 
LS 5039 appears as a bright - 11.2 star with an 07 V((f)) 
spectral type. Based on this evidence, M97 proposed the sys- 
tem to be a high mass X-ray binary (HMXRB). 

In addition, this system has been found to be active at radio 
wavelengths. Its radio counterpart (NVSS J182614- 145054) 
is a bright, compact and moderately variable radio source 
in excellent sub-arcsecond agreement with the optical star 

Send offprint requests to: M. Ribo, mribo@mizar.am.ub.es 



(Marti et al. 1998). All these facts point to the peculiar na- 
ture of RX J1826.2-1450/LS 5039, and suggest a classifica- 
tion among the selected group of radio emitting HMXRB. 

In order to explore how this source behaves compared to 
other members of its class (e.g. Cygnus X-1, LS Ih-61°303 and 
SS 433), we have analyzed the corresponding X-ray data from 
the All Sky Monitor (ASM) and the Proportional Counter Ar- 
ray (PCA) on board the satellite Rossi X-ray Timing Explorer 
(RXTE). In Sect. ^ we present an X-ray timing analysis based 
on both the ASM and the PCA instruments. The ASM data are 
suitable to study the long-term (days to months) temporal be- 
havior of the source, whereas X-ray variability on shorter time 
scales (seconds to hours) is better investigated with the PCA. 
In Sect. H a PCA spectroscopic analysis is presented, with the 
different spectral models that fit the data being examined and 
discussed. 

In the radio domain, RX J1826.2-1450/LS 5039 was in- 
cluded at our request in the list of radio sources routinely mon- 
itored at the Green Bank Interferometer (GBI)[} At the time of 
writing, the radio light curves cover ^ 4 months of observa- 
tions. In Sect. Q we present the GBI radio data so far acquired 
with some discussion on the source variability and spectral in- 
dex properties. Finally, we conclude in Sect. ^ with a brief com- 
parative discussion of RX J1826.2— 1450/LS 5039 versus other 
radio loud HMXRB. 

Hereafter, we will refer to the source as RX J1826.2— 1450 
when discussing the X-rays. In the optical/radio context the 
LS 5039 designation will be preferred. 

2. X-ray timing analysis 

2.1. The ASM/RXTE data 

The ASM database analyzed in this paper spans for more than 
two and a half years (1996 February-1998 November) and con- 
tains nearly 800 daily flux measurements in the energy range 
1.5-12 keV. Each data point represents the one-day average of 
the fitted source fluxes from a number (typically 5-10) of in- 
dividual ASM dwells, of ^ 90 s each (see Levine et al. 1996 

' The Green Bank Interferometer is a facility of the USA National 
Science Foundation operated by NRAO in support of the NASA High 
Energy Astrophysics programs. 
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Fig. 1. ASM one-day average light curve of RX J1826.2- 1450 in the 1.5-12 keV energy band. 
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for more details). The one-day average light curve is shown 
in Fig. |l]. The big gap between Modified Julian Date (MJD) 
~ 50400 and ^ 50500 corresponds to the passage of the Sun 
close to the source during the first year of observations. This 
gap repeats the following two years (near MJD 50800 and 
51150), but it happens to be less severe and a few flux mea- 
surements were then possible. 

Most of time, the source is at the threshold of ASM de- 
tectability. Nevertheless, we have searched for possible period- 
icities in the range from 2 to 200 d. The methods employed 
were the Phase Dispersion Minimization (PDM) (Stellingwerf 
1978) and the CLEAN algorithm (Roberts et al. 1987). Our ap- 
proach here is essentially the same as in Paredes et al. (1997) 
when analyzing the periodic behavior in the X-ray light curve 
ofLS lH-61°303. 

After applying both the PDM and CLEAN methods to the 
ASM data, a period of ^ 52.7 d stands prominently. This pe- 
riodicity corresponds to the detection of some kind of active 
events that appear rather evident at first glance in Fig. ^ Nev- 
ertheless, a careful inspection of the data reveals a suspicious 
detail. All those active events take place when the data by dwell 
coverage is rather poor (less than 5 dwells per day), thus re- 
ducing the statistical significance of the corresponding one-day 
average. For some instrumental reason, the ASM coverage be- 
comes poorer than normal every ^ 53 d or so and, in the case 
of a weak X-ray source like RX J1826.2— 1450, this can affect 
somehow the period analysis. Therefore, the ^ 52.7 d period is 
very likely to be an instrumental artifact. Indeed, after remov- 
ing all daily points resulting from less than 5 dwells 20% 
of total), the timing analysis reveals no significant period in the 
range from 2 to 200 d. 



Table 1. Log of the PCA observation of RX J1826.2-1450 



Date 


Start time 


Stop time 


Luminosity" 




(TT) 


(TT) 


(erg s"i) 


08/02/98 


01:05:33 


04:18:14 


6.4 X 10^" 


08/02/98 


20:07:34 


00:57:14 


5.3 X 10^* 


16/02/98 


18:42:25 


20:35:14 


5.3 X 10^* 



' in the energy range 3-30 keV and for an assumed distance of 3. 1 kpc 



2.2. The PCA/RXTE data 

Additional observations were made with the PCA instrument 
on 1998 February 8 and 16. The total on-source integration 
time was 20 ks. The PCA is sensitive to X-rays in the energy 
range 2-60 keV and comprises five identical co-aligned gas- 
filled proportional counter units (PCUs), providing a total col- 
lecting area of ^ 6500 cm^, an energy resolution of < 18 % at 6 
keV and a maximum time resolution of 1/xs. Our analysis was 
carried out in the interval 3-30 keV since the PCU windows 
prevent the detection of photons below 2.5 keV, whereas 
above 30 keV the spectrum becomes background dominated. 

Good time intervals were defined by removing data taken at 
low Earth elevation angle (< 8°) and during times of high par- 
ticle background. An offset of only 0.02° between the source 
position and the pointing of the satellite was allowed, to ensure 
that any possible short stretch of slew data at the beginning 
and/or end of the observation was removed. Table [l| shows the 
journal of the PCA observations, while the light curve of the 
entire observation is presented in Fig. ^. 

The PCA observations were used to study the time variabil- 
ity on various time scales. Continuous stretches of clean data 



M. Ribo et al.: X-ray and radio observations of RX J1826.2-1450/LS 5039 



3 



30 



25 



20 



15 



20 



15 



10 

25 



20 



15 



10 



h*4 


1 




' 


1 5000 


1 0000 


1 5000 


,1.1 

1 , , 


1 1 . . 

ItV * *t 

1 , , , 


1 11 11 _ 

1 , , , , " 


70000 


75000 


80000 


1 , , 



1 . . . 




, , , , _ 



755000 



760000 

Time (s) 



765000 



Fig. 2. 3-30 keV light curve of RX J1826.2- 1450 covering the 
entire PCA observation. Time is JD 2450 853.047 and the bin 

size is 180 s 



were selected from the light curve of the entire observation. To 
reduce the variance of the noise powers, these intervals were 
divided up into segments of 8192 bins each, with a bin size of 
10 ms. Then the power density spectra for each segment were 
calculated and the results averaged together Fig. |^ shows the 
characteristic power spectrum in the frequency range 0.01-50 
Hz. The dashed line represents the 95% confidence detection 
limit (van der Klis 1989). As it can be seen no power exceeds 
this value; the distribution of powers is flat at a level of 2, con- 
sistent with Poissonian counting statistics. Following van der 
Klis (1989), we can set a 95% upper limit of 60% on the rms of 
a pulsed source signal in the range 0.01-50 Hz. This relatively 
high limit is a consequence of the faintness of the source. On 
longer time scales, longer intervals 3200 s) were considered 
but no evidence for pulsations was found either. 

Likewise, we folded the light curve onto a set of trial peri- 
ods with the FTOOLS software package (a technique very sim- 
ilar to the PDM) and looked for a peak in the versus period 
diagram. None of the peaks found were statistically significant 
enough. Thus, we conclude that no coherent periodicities were 
detected in the range ^ 0.02 to ^ 2000 s. 

The mean X-ray intensity in the energy range 3-30 keV 
shows a slight decreasing trend with 19.7±0.2 count s^^ at the 
beginning of the observation (upper panel of Fig. ^) compared 
to 16.6±0.1 count s^^ and 16.2±0.2 count s^^ for the middle 
and bottom panels of Fig. ^ respectively. The fractional rms of 
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Fig. 3. Characteristic power spectrum of RX J 1826.2— 1450. 
The dashed line represents the 95% confidence detection Umit 



the 3-30 keV light curve corresponding to the entire observa- 
tion is 9%. 

The fact that no X-ray pulsations have been found in 
RX J1826.2— 1450 is consistent with the proposed idea that ra- 
dio emission and X-ray pulsations from X-ray binaries seem 
to be statistically anti-correlated (Fender et al. 1997), i.e., no 
X-ray pulsar has ever shown significant radio emission. 

3. X-ray spectral analysis 

3.1. Spectral fitting 

Since the light curve of the entire observation does not show 
sharp features, i.e., there is no significant spectral change 
throughout the observation, we obtained one average PCA en- 
ergy spectrum from the complete observation (Fig. 

Acceptable fits of the X-ray continuum were obtained with 
an unabsorbed power-law model, giving a reduced x^=l • 14 for 
56 degrees of freedom (dof). A multi-colour disk model, as ex- 
pected from an optically thick accretion disk (Mitsuda et al. 
1984) plus a power-law gave a reduced Xy=l-1 1 for 53 dof. The 
best-fit results are given in Table |[ Bremsstrahlung and two 
blackbody component models did not fit the data. Although the 
addition of a blackbody component to the power-law formally 
produces an acceptable fit, the value of the blackbody normal- 
ization was very low, with the error bar close to zero. In fact, 
an F-test shows that the inclusion of a blackbody component is 
not significant. 

The most salient feature that appears in the spectrum of 
RX J 1826.2- 1450 is a strong iron fine at - 6.6 keV (Fig. |). A 
Gaussian fit to this feature gives a line centered at 6.62±0.04 
keV, with an equivalent width (EW) of 0.75 ±0.06 keV and 
a FWHM of 0.9±0.2 keV. The high EW{Fe) value indicates 
that a large amount of circumstellar matter is present in the sys- 
tem. Unfortunately, the PCA energy resolution prevents from 
distinguishing between a broad line or two narrow components. 

A hydrogen column density of ^ (2^3) ^ ^0^^ cm^^ is 
found from the fit. This value is, however, not very well con- 
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Table 2. Spectral fit results for the power-law model plus a 
Gaussian component for the iron line. Uncertainties are given 
90% confidence for one parameter of interest. The spectrum 
was fitted in the energy range 3-30 keV 



Parameters 


Fitted values 


Nh (xlO^' cm"'') 




a 


1.95±0.02 


EnneiFe) (keV) 


6.62±0.04 


EWunem (keV) 


0.75±0.06 


FWH Munc(Pe) (keV) 


0.9±0.2 


o-iinc(Fe) (keV) 


0.39±0.08 


Xr (dof) 


1.14 (56) 



strained. In fact, it is consistent with zero. The difficulty in con- 
straining the hydrogen column density from our X-ray data can 
be attributed to the fact that the interstellar gas mainly absorbs 
X-ray photons with energies lower than 2-3 keV, i.e., outside of 
the energy range considered here. Nevertheless, it agrees with 
the value obtained from optical observations. From the A 4430 
and A 6284 interstellar bands M97 found E{B -V) = 0.8±0.2 
Using the relation A^h = 5.3 x 10' 
& Schmitt 1995), we obtain that - (4±1) x 



■21 cm-2 E{B - V) (Predehl 



which is consistent, within the errors, with the X-ray observa- 
tion value. 

3.2. A black hole or a neutron star? 

As mentioned above, no pulsations are found in the X-ray flux 
of RX J1826.2-1450/LS 5039. At the time of writing this 
paper approximately 55% of the optically identified massive 
X-ray binary systems are pulsars, rising to 67% if suspected 
HMXRB are included. About 75% of the identified and sus- 
pected X-ray pulsars have spin periods below 100 seconds. The 
detection of pulsations would rule out a black hole companion. 

Interestingly, unlike typical X-ray pulsars, the energy spec- 
trum of RX J 1826.2— 1450 shows no exponential cut-off at 
high energy despite that significant emission is seen up to 30 
keV. The absence of both, pulsations and high energy cut- 
off may indicate that the compact companion is a black hole 
rather than a neutron star. The persistent non-thermal radio 
emission of RX J1826.2— 1450/LS 5039 is also reminiscent, 
among others, of the classical black hole candidate (BHC) 
Cygnus X- 1 . However, the data are not conclusive as counter- 
examples of these properties can be found. For example, the 
system LS 1h-61°303, which seems to contain a neutron star 
and does not exhibit pulsations or cut-off in the 10-30 keV 
range. The multi-colour disk model, although formally fitting 
the data, does not help either. First, given the low luminosity 
(Lx <10"^^ erg s^^ in the energy range 2-10 keV) the sys- 
tem would be in the so-called low-state. We would not ex- 
pect then to detect a strong soft component. Second, the fit 
provides an unrealistic value of the disk internal radius of 
i?in cosi/2(6') ~ 0.3 km. Finally, the lack of detected pulsa- 
tions may be just due to the faintness of the X-ray emission in 
view of the rather high upper limit found in Sect. 2.2 for pulsed 
emission. 

Unfortunately, the source is too faint at energies above 30 
keV to be detected with the HEXTE instrument. Thus, we can- 
not confirm from the present data whether the hard tail that 
characterizes the energy spectrum of BHCs at high energies is 
indeed present. In any case, the issue of a possible black hole 
in RX J1826.2-1450/LS 5039 is likely to be set in the future 
by obtaining a spectroscopic mass function in the optical. 

4. Radio observations 

4.1. Radio variability properties ofLS 5039 

The radio observations reported here consist of daily flux den- 
sity measurements with the GBl within the GBl-NASA Moni- 
toring Program. The source was observed at the frequencies of 
2.25 and 8.3 GHz. In Fig. ||, we show the corresponding flux 
density light curves during the first ^ 4 months of monitoring 
so far available. The bottom panel displays the spectral index 
a (where Si, oc v"^) computed between these two frequencies. 

The source was found to be detectable at radio wavelengths 
throughout all the time. This behavior is better evident at 2.25 
GHz where the source is brighter. The average GBl flux densi- 
ties and their respective rms are 5*2. 25 ghz — 31(±5) mJy and 
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Ss.3 GHz = 14(±5) mJy. The typical day-to-day variability in 
the GBI data does not exceed ^ 30 %. There may be however 
some exceptions, such as for example around Modified Julian 
Date (MJD) 51075, 51086 and 51176. Here, the flux density of 
LS 5039 seems to have varied by more than a factor of ~ 2 on 
less than one day. Excluding these episodes, the source never 
exhibited well defined radio outburst events. In general terms, 
the observed radio behavior confirms the early suggestions by 
Mart! et al. (1998) concerning the persistent and moderately 
variable nature of the radio emission. 

Timing analysis of radio light curves has proven to be in 
some cases a useful tool to detect orbital periods. For exam- 
ple, the orbital period signatures of 26.5 and 5.6 d are visible in 
LS Ih-61°303 and Cygnus X-1, respectively (Taylor & Gregory 
1984; Pooley et al. 1998). We have thus searched for long-term 
periodicities in the GBI data of LS 5039. Given the span of the 
radio observations, the search was restricted between 2 and 50 
d. The methods used were again the PDM and CLEAN, men- 
tioned in Sect. 



2.1 



Unfortunately, no convincing period was 
detected in this process. A longer time span is likely to be re- 
quired before a reliable search can be attempted for this rela- 
tively weak radio source (specially at 8.3 GHz). 



4.2. Non-thermal radio spectrum and brightness temperature 

From the GBI data, the weighted average spectral index is 
found to be a = — O.S^o's- This value is also in good agree- 
ment with the results by Marti et al. (1998) obtained a few 
months before, thus suggesting that the non-thermal radio spec- 
trum is a persistent property of the source. 

In addition to negative spectral indices, the brightness tem- 
perature estimates for LS 5039 clearly yield to non-thermal val- 



ues, hence supporting the mechanism of synchrotron radiation 
for that source. The apparent one-day variability observed in 
the GBI data around MJD 51075, 51086 and 51 176 would im- 
ply, from light time travel arguments, that the emitting region 
is smaller than about 2.6 x 10^^ cm. If we assume a 3.1 kpc 
distance to the source, the corresponding angular size is found 
to be 6* < 0"05, yielding a lower limit of Th > 6 x 10^ K (at 
2.25 GHz). This lower limit is not far from the Tf, > 4 x 10^ K 
estimate by Marti et al. (1998), based on the unresolved nature 
of the source with the VLA. 



5. Comparison to other radio loud HMXRB 

SS 433, LS lH-61°303 and the BHC Cygnus X-1 are classical 
examples of HMXRB with detectable radio emission (Penninx 
1989), and all of them are also under GBI monitoring. In or- 
der to facilitate the comparison of LS 5039 to these sources, 
we have summarized some relevant parameters in Table |[ 
They include: the quiescent X-ray luminosity, as derived from 
ASM/RXTE data (extrapolated from PCA/RXTE in the case of 
LS 5039 because the signal to noise ratio is much higher); the 
weighted average radio luminosity and spectral index in qui- 
escence, both based on GBI data; the nature of the compact 
object; the spectral type of the companion; the orbital period of 
the binary system and the distance to the source. 

As it can be seen in Table|^, the radio luminosity of LS 5039 
is very similar to that of Cygnus X-1. This is, of course, pro- 
vided that the distance adopted is correct. Their respective 
spectral indices seem to be intrinsically different, both sources 
being persistent at radio wavelengths. Cygnus X-1 also ex- 
periences strong X-ray variability due to changes in its state, 
that we are not aware of in our source. The LS Ih-61°303 ra- 
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Table 3. Average properties of radio emitting HMXRBs from ASM/RXTE and GBI data 



HMXRB 


Lx(1.5-12keV) 


irad(0. 1-100 GHz) 


Q 


Compact object 


Companion 


Orbital 


Distance 




(erg s"^) 


(erg s"^) 


(S. oc I/") 






period (d) 


(kpc) 


LS 5039 


~ 5 X 10^" 


1.0 X 10^^ 


-0.5 


? 


07V((f)) 


? 


3.1" 


Cygnus X- 1 


~ 8 X 10^*^ (*) 


1.1 X 10^^ 


0.1 


black hole 


09.7 lab 


5.6 


2.5'' 


LS 1+61° 303 


~ 4 X lO^"* 


0.9 X 10^^ (*) 


-0.3 (*) 


neutron star ? 


BOVe 


26.5 


2.0^= 


SS 433 


~ 7 X 10^^ (*) 


3.2 X 10^^ (*) 


-0.7 (*) 


neutron star ? 


OB? 


13.1 


4.8'^ 



(*) Only data during quiescence has been considered 
" Motch et al. 1997 
^ Penninx 1989 
" Frail & Hjellming 1991 
Vermeulen et al. 1993 



die properties during quiescence are also very comparable to 
those of LS 5039, as well as their respective X-ray luminosi- 
ties. In contrast, LS I+61°303 undergoes radio outbursts every 
26.5 d, the orbital period of the system, while LS 5039 never 
had strong outbursts during the GBI observations. The X-ray 
and radio luminosities of SS 433 in quiescence are much higher 
than those of LS 5039. However, we notice that their Lrad/^x 
ratios are practically the same. 

In general terms, the average properties of LS 5039 do not 
deviate extraordinarily from those of other radio loud HMXRB. 
Since even the well accepted members of this class are not an 
homogeneous group, the belonging of LS 5039 to this category 
appears as very plausible to us. 

6. Conclusions 

We have presented a general overview of the X-ray and 
radio emission properties of the massive X-ray binary 
RX J 1826.2- 1450/LS 5039. Our X-ray and radio results are 
mostly based on long term (few months) monitorings of the 
source, with our main conclusions being: 

1. In the X-rays, a timing analysis has been performed show- 
ing neither pulsed nor periodic emission on time scales of 
0.02-2000 s and 2-200 d, respectively. The X-ray spec- 
trum has been found to be significantly hard (up to 30 keV), 
with no cut-off required. It can be fitted satisfactorily with 
a power-law plus a strong Gaussian iron line. 

2. At radio wavelengths, the GBI monitoring confirms the 
long-term persistence of the RX J1826.2- 1450/LS 5039 
radio emission in time scales of months, always with a non- 
thermal synchrotron spectrum. The day-to-day variabiUty 
continues to be moderate most of the time (^ 30 %), and 
no strong radio outbursts have been observed. 

3. The classification of RX Jl 826.2- 1450/LS 5039 among 
the radio loud HMXRB group is reinforced. Although some 
specific differences with other members of this class do ex- 
ist, noticeable similarities can be found. 
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